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ABSTRACT 


Simulated surface conditions of the Goddard Earth Observing System model, version 5 
(GEOS-5) atmospheric general circulation model (AGCM) are examined for the contemporary 
Greenland Ice Sheet (GrIS). A surface parameterization that explicitly models surface processes 
including snow compaction, meltwater percolation and refreezing, and surface albedo is found to 
remedy an erroneous deficit in the annual net surface energy flux and provide an adequate 
representation of surface mass balance (SMB) in an evaluation using simulations at two spatial 
resolutions. The simulated 1980-2008 GrIS SMB average is 24.7±4.5 cm yr -1 water-equivalent 
(w.e.) at V 2 0 model grid spacing, and 18.2±3.3 cm yr -1 w.e. for 2° grid spacing. The spatial 
variability and seasonal cycle of the V 2 0 simulation compare favorably to recent studies using 
regional climate models, while results from 2° integrations reproduce the primary features of the 
SMB field. In comparison to historical glaciological observations, the coarser resolution model 
overestimates accumulation in the southern areas of the GrIS, while the overall SMB is 
underestimated. These changes relate to the sensitivity of accumulation and melt to the resolution 
of topography. The GEOS-5 SMB fields contrast with available corresponding atmospheric 
models simulations from the Coupled Model Intercomparison Project (CMIP5). It is found that 
only a few of the CMIP5 AGCMs examined provide significant summertime runoff, a dominant 
feature of the GrIS seasonal cycle. This is a condition that will need to be remedied if potential 
contributions to future eustatic change from polar ice sheets are to be examined with GCMs. 
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1. Introduction 


The Greenland Ice Sheet (GrIS) influences global climate in part through the storage and 
release of freshwater into the ocean. Changes in mass of the grounded ice sheet resulting in 
freshening of the adjacent ocean through runoff and discharge have implications for global sea 
level and for North Atlantic thennohaline circulation (e.g., Fichefet et ah, 2003; Gerdes et ah, 
2006). In recent years, the GrIS has made a positive contribution to global sea level rise (e.g., 
Hanna et al., 2013; Bindoff et al., 2007; Shepherd et al., 2012). Vaughan et al. (2013) 
characterize as “very likely” that the rate of ice loss from the GrIS over the last two decades has 
accelerated, and that the loss is equally partitioned between outlet glacier discharge and surface 
melt (van den Broeke et al., 2009). Surface mass balance (SMB)- the balance of accumulation 
minus surface wastage terms- is thus an important variable for the evaluation of GrIS conditions. 
Along with surface temperature, SMB is a fundamental input field for dynamical ice sheet 
models (ISMs) used in prognostic assessments of GrIS total ice volume and glacier discharge 
(Bindschadler et al., 2013; Nowicki et al., 2013; Alley and Joughin, 2012). In the evaluation of 
ISMs, uncertainty associated with external surface forcings has been found to be as important as 
uncertainty produced by the different fonnulations of the models themselves (Nowicki et al., 
2013). 

Interest in examining the surface representation of the GrIS in general circulation models 
(GCMs) is motivated by a need to assess changes to its eustatic contribution under potential 
climate scenarios and to understand how the partitioning of the contribution might change with 
time (e.g., Parizek and Alley, 2004). Skill in reproducing GrIS SMB is known to vary widely 
among GCMs. For example, Smith (1999) found that “realistic” estimates of SMB terms were 
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obtainable from integrations with a l 7 /* 0 grid interval (approximately 200km x 70km), while 
Murphy et al. (2002) and Bengtsson et al. (2011) found that simulated GrIS SMB became 
unrealistically negative for grid spacings of 125km or greater. Two issues of relevance for GCM 
simulations are (1.) the spatial resolution afforded by the model grid and (2.) the extent to which 
physical processes are represented. With regards to spatial resolution, Meehl et al. (2007) 
observed that GCMs surveyed in the fourth assessment report of the Intergovernmental Panel on 
Climate Change (IPCC AR4) were unable to produce a “good simulation” of SMB, in part 
because the available spatial resolution was insufficient to adequately resolve the coastal 
escarpment of ice sheets and its associated orographic effects on precipitation and melt (e.g., 
Glover et al., 1999). In the fifth assessment report (AR5), Flato et al. (2013) noted that 
agreement with observations is limited for the polar ice sheets. Church et al. (2013) further 
indicated that resolution of a few tens km or finer is necessary to resolve strong gradients in the 
SMB field. Efforts have been made to apply downscaling approaches for producing a higher 
resolution SMB field from GCM models (Hanna et al., 2005; Agosta et al., 2013; Bougamont et 
al., 2005; Jarosch et al., 2012; Mernild and Liston, 2012). These strategies are advantageous for 
coupling with ISMs, which have spatial resolutions of a few km or less (Pollard et al., 2010; 
Lipscomb et al., 2013; Ridley et al., 2005; Vizcaino et al.; 2010). 

With regards to physical processes including surface albedo, snow ageing, and the 
treatment of melt, IPCC reports indicate deficiencies in the representation of these phenomena in 
GCMs. Deficiencies in SMB fields produced by AR4 models were found to be a result of 
limitations in the land surface schemes employed, which lacked a representation of melt water 
refreezing within the snow pack and snow albedo variations (Meehl et al., 2007). In addition, the 
effects of GrIS meltwater runoff on the North Atlantic meridional overturning circulation were 
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commonly not included in AR4 simulations (Randall et al., 2007). More recently, the AR5 report 
also noted that ice sheet processes were not generally well represented in climate models 
(Bindoff et ah, 2013). Several recent studies have advanced more sophisticated surface schemes 
as an essential step for an improved SMB representation (Punge et al.,2012; Vizcaino et al., 
2013). As surface processes on the GrIS have characteristic length scales and occur in particular 
regions (Wild et al., 2003; Quiquet et al., 2012), spatial resolution and model complexity are 
significantly linked. This needs to be taken into consideration for an evaluation of the ice sheet 
surface representation. 

In this study, an improved surface representation is evaluated in atmospheric GCM 
(AGCM) integrations of the contemporary climate at both coarse and high spatial resolution. The 
Goddard Earth Observing System model, version 5 (GEOS-5) is assessed in comparison to a 
standard model configuration, available observational datasets, and other GCMs. The paper is 
organized as follows. Section 2 describes the GEOS-5 model and the implementation of Stieglitz 
snow hydrology model over glacial surfaces. Datasets used for assessing model performance are 
also described. An appraisal of the effects of the snow model on the general climate in 
comparison to the control simulation are provided in section 3, as well as a detailed assessment 
of the resulting surface energy balance and SMB. A discussion of the results is given in 
section 4. 

2. Model Description, Datasets, and Method 

GEOS-5 is an AGCM maintained by the National Aeronautics and Space Administration 
Global Modeling and Assimilation Office (NASA GMAO; Rienecker et al., 2008; Molod et al., 
2012). The GEOS-5 AGCM is used for decadal climate prediction studies and as the atmospheric 
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component in coupled model simulations (Ham et al., 20 12a, b), routine numerical weather 
prediction (Reale et al., 2009), and serves as the background atmospheric model in the data 
assimilation system used for the Modern-Era Retrospective Analysis for Research and 
Applications (MERRA; Rienecker et al., 2011). The model employs a finite-volume dynamical 
core (Lin, 2004) that is integrated with various physics packages. In the version of GEOS-5 used 
for this study, a 30 arc-second version of the global elevation dataset of the Shuttle Radar 
Topography Mission (SRTM30; Farr et al., 2007) is supplemented by data over polar ice sheets 
(Bamber et al., 2001; Liu et al., 2001). 

The GEOS-5 land surface model is catchment-based, wherein a topographical index 
denotes a set of tiles occupying fractional areas below each atmospheric grid column (Koster et 
al., 2000). These tiles are used to indicate land surface types that represent sub-grid scale 
heterogeneity in prognostic hydrologic variables, such as soil moisture. The index of each sub- 
grid scale catchment tile is determined from a high resolution raster- currently from the 5km 
resolution dataset of the USGS Global Land Cover Characterization project (Loveland et al., 
2000). Fluxes computed on surface tiles are related to the atmospheric column via an exchange 
grid. 

The terrestrial (i.e., non-glacial) land surface model in GEOS-5 employs a modified 
version of the Stieglitz snow hydrology model, which provides an explicit representation of snow 
densification, meltwater runoff, percolation, and refreezing, and surface albedo (Lynch-Stieglitz, 
1994; Stieglitz et al., 2001). In the adaptation for GEOS-5, the Stieglitz model has three layers 
when snow cover is present, with a maximum thickness of the top layer not exceeding 8cm. For 
each layer, prognostic variables of heat content, equivalent snow water content, and snow depth 
are maintained. A volumetric water holding capacity is denoted for each layer. Snow compaction 
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— 3 

is based on the parameterization of Kojima (1967). A maximum snow density of 500 kg m is 
prescribed, with snow mass exceeding this value becoming runoff. The snow model includes 
representations of processes including snowpack growth and ablation, melting, refreezing, and 
sublimation. For GEOS-5, snow albedo is prescribed based on a linear relation with snow 
density in the uppermost layer, with reductions imposed by vegetation masking and fractional 
snow cover. For a given model time step, the land surface scheme solves for the energy balance 
at both the snow-atmosphere and the land-snow interfaces. 

On grounded ice, hydrologic processes including runoff were not considered in the 

original model configuration; hence, an ice-sheet SMB was not produced. Ice sheets were 

represented as a 7 cm water-equivalent (w.e.) surface ice layer (Fig. 1). The temperature 

corresponding to this surface ice layer was determined from turbulent, radiative, and subsurface 

energy fluxes using semi-implicit time stepping. In the original model configuration, the 

subsurface energy flux was determined from the prognostic surface layer temperature and a fixed 

temperature of 230 K at 2m depth, and the land-ice surface albedo was fixed at 0.775. These 

values were applied in early model development, but are thought to have been selected as 

representative of the climate of polar ice sheets. In the semi-implicit scheme, atmosphere and ice 

surface temperatures respond mutually to these parameterizations. Cullather and Bosilovich 

(2012) examined the surface energy budget in MERRA over Greenland and Antarctica. The 

subsurface energy parameterization was found to produce spurious annual mean flux imbalances 
_2 

of up to 30 W m “ locally. The motivation for addressing limitations in the surface representation 
over grounded ice in GEOS-5 thus included the need to remedy existing limitations in the 
surface energy budget representation as well as for producing a realistic estimate of runoff for 
use in a coupled atmosphere-ocean-dynamical ice sheet model configuration. 
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Over glaciated surfaces, snowfall typically exceeds summer melt. Snow hydrology 
representations commonly used in GCMs over non-glaciated surfaces accumulate snow without 
a mechanism for converting it to glacial ice and, over time, would produce a snowpack over ice 
sheets with nonphysical depth. While this does not constitute a direct impediment to the 
perfonnance of the atmospheric model, the conducted heat flux and other properties may not be 
adequately parameterized or resolved with available resources through an unrealistically deep 
snow pack, thus having an adverse effect on simulated surface energy fluxes. To compensate, the 
use of a snow hydrology model over ice sheets is often prohibited as in the original version of 
GEOS-5, or a cap on the size of the snowpack is provided (e.g., Lipscomb et ah, 2013). Snow 
mass exceeding the cap in snow depth or in density is then directly routed to the oceans. This 
excess term is often referred to as “frozen runoff’, and may be seen as a poor man’s simulation 
of ice calving. The criteria necessary for an improved land ice surface representation may then be 
seen as one that maintains consistency with the snow hydrology representation used on other 
terrestrial surfaces while accommodating specific issues associated with snowfall on ice sheets. 

The new configuration for land ice surfaces (Fig. 1(b)) indicates an implementation of the 
snow hydrology model with a cap on the maximum snow depth at 15 m. This depth is 
comparable with other models (e.g., Bugnion and Stone, 2002). In practice, snow depth is limited 
by the Stieglitz model maximum value for snow density and is typically less than half the 
imposed maximum depth at any location. In this elemental approach, the Stieglitz model 
operates as it does over non-glaciated surfaces, but with an increased vertical resolution of 
15 layers to adequately simulate the heat flux through the snowpack. The upper 5 layers of the 
snowpack are constrained in depth to resolve the near-surface temperature gradient, with the top- 
most snow layer restricted to a maximum thickness of 8 cm. Snow cover over ice is allowed to 


8 


158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 


be fractional, and melting and runoff may occur on exposed bare ice surfaces. A 15-layer ice 
column is applied for the conduction of heat below the snow-ice interface, with a zero heat flux 
condition imposed at the lower boundary. Surface radiative properties for the snowpack are taken 
from Greuell and Konzelmann (1994), in which albedo is a linear function of snow density in the 
uppermost snow layer. In this relation, albedo decreases more rapidly with increasing snow 
density than is prescribed by Stieglitz et al. (2001). This is to account for the higher density of 
fresh snow found over ice sheets (Greuell and Konzelmann, 1994). Recent studies have a found 
more realistic representation of snow albedo using parameterizations of snow grain size rather 
than snow density as used in GEOS-5 (van Angelen et al., 2012). In particular, a density- 
dependent albedo scheme was found to underestimate albedo in wet snow conditions. 
Nevertheless, credible results are obtainable with snow density parameterizations (e.g., Ettema et 
al., 2009; Ettema et al., 20 10a, b; Slater et al., 1998). A value of 0.6 is used for bare ice albedo. 
Additional characteristics of the snow model are as previously described in Lynch-Stieglitz 
(1994) and Stieglitz et al. (2001). Thus for the new configuration, snow hydrology over land ice 
remains consistent with that used over non-glaciated surfaces, with the exception of increased 
vertical resolution for larger snow depths and modification of the snow albedo parameterization. 

The approach of this study is to evaluate the performance of the GrIS surface 
representation in AMIP-style simulations of GEOS-5 forced with sea ice and SST fields for the 
period 1979-2008 (AMIP: Atmospheric Model Intercomparison Project; Gates et al., 1998; 

SST: sea surface temperature; Reynolds et al., 2002). Two simulations are run at a horizontal 
resolution of 2° latitude by 2 l A° longitude and 72 hybrid-sigma coordinate vertical levels. Two 
simulations of V 2 0 x S A° are also presented and will be the main focus for comparisons with 
in situ observations. Hereafter, these simulations are referred to as GSN2 for the 2° integrations 
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using the snow hydrology model over land ice regions, GSNV 2 for the V%° run using the snow 
hydrology model, and CNTRL2 and CNTRLV 2 for the control simulations at 2° and Vi°, 
respectively. Figure 2 shows the differences in the spatial representation of Greenland 
topography. While model GrIS elevations are both initially based on Bamber et al. (2001), the 
process leading to model topography that is appropriately scaled for each GCM spatial resolution 
is complex. The V 2 0 model topography closely approximates that of Bamber et al. (2001). The 2° 
model topography is lower over the central and western areas of the GrIS, but the 50m contour 
encompasses a larger area. Surface characteristics in these simulations are compared with the 
control simulation, with contemporary GCM and regional model output, and with available 
observations and gridded datasets. These include available in situ observations of surface 
temperature, energy fluxes, and glaciologically-derived SMB, and with gridded fields of 
satellite-derived temperature and prognostic SMB from regional climate models and from the 
European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis Interim product 
(ERA-I: Dee et al., 2011). Monthly fields of the ERA-I were obtained for the years 1980-2008 at 
a regular grid spacing of 0.75° x 0.75°. A collection of 15 automatic weather stations (AWS) 
known as GC-Net are distributed on the GrIS and report hourly measurements of temperature, 
wind, surface radiation fluxes, and other variables (Steffen and Box, 2001). Twelve AWS 
stations from the Danish Meteorological Institute (DMI; Carstensen and Jorgensen, 2011) 
located on non-glaciated land surfaces on the GrIS periphery are also used. Additional datasets 
are described in context below. 

3. Results 

a. Surface energy budget 
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From the perspective of the atmospheric model, the net surface energy flux is the sum of 

radiative flux components and turbulent fluxes of sensible and latent heat. In the absence of melt 

processes, it is in agreement with the ground flux through the ice. For the control model version 

of GEOS-5, the annual mean simulated net surface energy flux over the GrIS consists of a large 

downward (negative) flux of greater than 25 W m "in coastal regions and greater than 5 W m 

over the ice sheet interior. This result from AMIP simulations is similar to MERRA (Cullather 

and Bosilovich, 2012). The flux results from the fixed sub-surface temperature parameterization 

of 230K. While the parameterization was thought to provide a more realistic diurnal cycle, the 

annual mean surface flux magnitudes are erroneous. Net surface energy fluxes over ice sheets 

may result from geothermal sources but are thought to be small (e.g., Reijmer and Oerlemans, 

2002) and are not explicitly simulated here. The annual net surface flux may also differ from 

zero in locations of surface melt or refreezing. For the new scheme depicted in Fig. 1(b), the net 

surface flux is near zero for the bulk of the GrIS with the exception of small coastal regions. In 

_2 

particular, magnitudes greater than (-) 10 W m are found only on the southwestern edge of the 
GrIS which are associated with surface melt. For the new scheme, a net heat flux of zero is not 
guaranteed over non-melting land ice surfaces, due to the removal of mass associated with 
capping the snow density and depth, but this heat flux due to “frozen runoff’ is small. 

Hoch (2005) obtained radiative flux values at Summit (73°N, 38°W) for the period July 
2000 to July 2002, and turbulent flux values from a 50m tower for the period June 2001 to July 
2002. Components of the surface radiative flux were also obtained from AWS stations as part of 
the Programme for Monitoring of the Greenland Ice Sheet for the period 2008-2010 (PROMICE; 
van As and Fausto, 2011), and as part of the Greenland Analogue Project along the 
Kangerlussuaq transect (K-transect) on the western GrIS near 67°N for the period 2003-2010 
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(van den Broeke et al., 2008). K-transect site 9 (S9) was located within 91 km of the western ice 
margin. The PROMICE AWS stations are also typically located along the periphery of the GrIS. 
Comparisons of the V 2 0 model with Summit observations shown in Fig. 3(a) indicate 

_2 

disagreements with the surface net flux from the CNTRLV 2 simulation of greater than 20 W m 

_2 

in summer months, and an annual mean difference of 8 W m . This summertime bias in the net 

flux is eliminated with the new surface representation. The improvement results from a better 

treatment of the shortwave flux, and error cancellation. As seen in Fig. 3(b), differences with the 

July observed net shortwave flux are reduced, from 20 VV in for the CNTRLA simulation to 

1 1 W m in the GSNA run. An important component of the improved representation is the 

Greuell and Konzelmann (1994) surface albedo parameterization. The application of the surface 

model allows the surface albedo to develop spatial heterogeneity, as shown in Fig. 4 for July. A 

considerable area of the ice sheet above the ablation zone has simulated albedos greater than the 

0.775 prescribed value in the control simulation. The spatial distribution shown in Fig. 4 is 

comparable to that given by Vizcaino et al. (2013) for RACM02 and the Community Earth 

System Model (CESM). In contrast to improvements found in the shortwave flux representation, 

_2 

an underestimate of the annual mean net longwave flux of 12 W m is increased in the GSNV 2 
_2 

run to 16 W m . The turbulent flux values from GSNA are not significantly different from the 
control. 

For K-transect station S9 (Fig. 3(c) ) located within the ablation and percolation zone, a 

_2 . . . _2 

net surface flux bias of 18 W m in CNTRFV 2 winter is reduced to 2 W m in the GSN!4 run. 
As with Summit, there is better agreement between the observed and new model simulated net 
shortwave flux in summer and autumn as compared to the control (Fig. 3(d) ). But improvement 
in the net flux in winter largely results from changes in turbulent fluxes as compared with the 
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_2 

control simulation. A sensible heat flux bias of 32 W m in January in the control model is 

reduced to 20 W m for the GSN'A model. However differences between model and observation 

for the net longwave flux are again larger with the new surface parameterization. The annual 

mean bias in the net longwave flux is 2 W m 2 for CNTRLA but 15 W m 2 for GSN'A 

Comparisons with other K-transect and PROMICE AWS locations indicate similar biases in the 

net longwave flux. In comparison to four PROMICE stations located on the GrIS, there is an 

”2 

underestimate of the modeled annual net longwave flux for the GSNV 2 model of 3 to 9 W m . 
These biases are found to result from an underestimate of the downwelling component. For the 
GSN'/i simulation, surface temperatures are improved in comparison to the CNTRLA model 
such that a compensating bias in the upwelling longwave flux is reduced, and an underestimate 
in the downwelling flux remains. This is related to the underestimation of clouds, atmospheric 
humidity, and/or transmissivity, and is common to most atmospheric models. 

Figure 5 presents a comparison of averaged near-surface air temperatures from GC-Net 
and DMI stations with model values from the V 2 0 simulations. For winter months, temperatures 
from the CNTRLA and GSNA simulations compare similarly, and the skill is comparable to 
other climate models (e.g., Rae et ah, 2011). In both simulations, near-surface DJF air 
temperatures correlate well with GC-Net (r = 0.98) and coastal DMI stations (r = 0.96). The 
model is biased by -2.6°C in comparison to GC-Net and +1.1 °C in comparison to DMI. 

Important differences occur between the two simulations in summer months. For the 
interior GC-Net stations, the average JJA temperature bias is -4.4°C for the CNTRLA 
simulation, but -1.4°C for the GSN'A run. Farger differences between the simulations occur for 
locations with warmer observed temperatures. For locations warmer than -5°C, the GSNA 
simulation is 3.8°C warmer than the control. For locations observed colder than -10°C, the new 
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scheme is 2.8°C warmer. This implicates the use of a cold, fixed subsurface temperature in the 
control simulation. The correlation between summer observed and simulated values is also 
slightly improved with the new parameterization, from r= 0.96 for the control to r = 0.97 for the 
new parameterization. For DMI stations, the average bias is reduced from -2.6°C in the control 
to -1.7°C with the new scheme, and the correlation has increased from r= 0.80 to r = 0.84. 

The spatial agreement may be further evaluated using monthly temperature data derived 
from the Moderate Resolution Imaging Spectroradiometer (MODIS; Hall et ah, 2012), which 
were obtained for the period 2003-2010. For comparison with MODIS, the model temperature 
was computed from the clear-sky emitted longwave flux field from the surface. The difference 
between the summertime model-computed temperature and MODIS for the CNTRL14 simulation 
is shown in Fig. 6(a). The figure shows magnitudes greater than (-) 2°C over most of the GrIS, 
with largest differences of greater than (-) 4°C inland from the southeastern escarpment and 
greater than (-) 5°C adjacent to the Scoresby Sund. Some edge effects are also apparent, 
particularly in the northern regions. With the new surface representation shown in Fig. 6(b), 
these differences are greatly reduced over much of the GrIS, with values colder than MODIS by 
more than 2°C in the southeast, and more than 3°C warmer along the western edge of the ice 
sheet. MODIS ice surface data are kn own to have a cold bias in comparison with thermochron- 
derived surface temperatures (e.g., Hall et ah, 2012; Koenig and Hall, 2010). Nevertheless the 
comparison shown in Fig. 6 indicates a substantial correction to surface temperatures with the 
new surface scheme. 

The evaluation of subsurface conditions is even more challenging, owing to a lack of 
reliable contemporary observational data. Observations from the GC-Net AWS network include 
subsurface temperature profiles produced from a thermocouple string. Such measurements are 
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problematic, as the vertical position of the string may change with time in an unknown manner 
from accumulation or settling (Steffen and Box, 2001). Nevertheless, a comparison of the 
modeled annual temperature wave with observations at two representative stations (Fig. 7) show 
similarities with the observed rapid damping with depth. For a low accumulation site, the 
simulated 10m depth annual cycle magnitude at TUNU-N (78° 01'N, 33° 59’W) of 0.7°C in 
GSN'A is comparable with the observed value of 0.5°C. This location may be contrasted with 
Swiss Camp (69° 34’N, 49°20'W), where interannual variability is large (Steffen and Box, 2001). 
Here the 10m depth annual temperature wave magnitude in GSNA of 1.2°C compares with 0.7°C 
from observations, and there is a bias of 2.5°C in the mean annual temperature. As well as 
measurement uncertainties, differences in the simulated magnitude of the annual cycle in the 
upper 3 to 4m may also result from comparison of point observations with the AGCM grid box. 
Given the uncertainties, such differences are not unexpected (e.g., Bugnion and Stone, 2002). 

The differences in the surface energy budget resulting from the new surface 
representation affect the local and regional atmospheric circulation. Shown in Fig. 8(a) are the 
pan-Arctic temperature differences between the CNTRLA and GSNA simulations for summer 
months. Note that again the largest temperature differences on glaciated land surfaces occur in 
locations where the surface values differ significantly from the 230K parameterized ice value of 
the control simulation. These locations are along the periphery of the GrIS. Smaller glaciated 
surfaces are apparent, including Ellesmere, Svalbard, and the Juneau Icefield. The difference in 
the pan-Arctic sea level pressure field shown in Fig. 8(b) indicates a local deepening of the 
pressure field in the North Atlantic near Greenland of about 1.5 hPa. The differences of the two 
model simulations are significant but are less than the observed interannual variability at coastal 
stations in southeastern Greenland. For example, the standard deviation of the observed JJA sea 
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level pressure at Tasiilaq (65.6°N, 37.6°W) is 2.3 hPa for the period 1980-2008. The pressure 
differences reflect a dynamical response to wanning of the ice sheet periphery. These differences 
also appear to influence the wave-2 pattern of the Northern Hemisphere, with a significant 
lowering of the pressure field over Siberia. 

b. Surface mass balance 

The ice-sheet SMB may be approximated by the net of three terms: precipitation- the 
vertical flux of solid and liquid water phases between the atmosphere and the surface; 
evaporation- the net vertical flux of water vapor between the surface and the atmosphere; and 
runoff- the net horizontal divergence of liquid water at the surface. Other terms including 
blowing snow are considered negligible here (e.g., Loewe, 1970) although studies have estimated 
substantial values in winter for the GrIS (e.g., Lenaerts et ah, 2012) and Antarctica (e.g., Budd et 
ah, 1966). 

Figure 9 shows the averaged annual SMB from GEOS-5 using the new surface 
representation at two different spatial resolutions. Again, the control version of GEOS-5 does not 
produce SMB. In the land surface scheme used in GEOS-5, surface characteristics for each 
fractional land surface type are aggregated to form a flux for the atmospheric model. For 
comparison with other models, the GEOS-5 fields shown in Fig. 9 are produced from the 
assembled global field, which includes non-glaciated land surfaces (runoff over oceans is zero). 

Also shown for comparison are the regional climate model simulation of RACM02 
1 1 km integration (Ettema et ah, 2009) interpolated to 5km grid spacing (Bindschadler et ah, 
2013), and the average of accumulated 12-hr prognostic integrations from ERA-I. The ERA-I 
surface fluxes denote the average of short-tenn model output initialized with assimilated 
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observations, and have been shown to better agree with in situ observations than other analyses 
(e.g., Chen et ah, 2011). The ECMWF model represents snow cover as one layer with a 
maximum depth of 7 cm, and meltwater percolation and refreezing processes are not represented 
(ECMWF, 2007; Ettema et ah, 2010b). The ERA-I prognostic snow albedo follows the 
parameterization of Douville et al. (1995) and is a function of snow age (ECMWF, 2007). ERA-I 
analyses are also used as lateral boundary conditions for the RACM02 regional model. 

RACM02 simulations are considered a state-of-the-art representation of SMB due to the use of 
high spatial resolution and representation of physical processes, and have been evaluated against 
in situ and satellite observations (e.g., Hanna et al., 2011; Shepherd et al., 2012). Results from 
the Modele Atmospherique Regional (MAR; Franco et al., 2012) are also cited. Church et al. 
(2013) characterized regional climate models as “the primary source” of ice-sheet SMB 
projections. The ERA-I SMB shown in Fig. 9 at 0.7° grid spacing presents the basic large-scale 
features for the GrIS as they are known, but with differences in the distribution of ablation zones 
as compared to RACM02. Along the western escarpment of the GrIS, local maxima near 
Qaanaaq and the Uummannaq Fjord system seen in the RACM02 field are not present in the 
ERA-I. The topographies of the given models and reanalyses are also indicated in Fig. 9. For 
RACM02, the Bamber et al. (2001) topography is contoured. 

The GSNA simulation closely approximates the RACM02 field in showing largest 
values of greater than 200 cm yr -1 w.e. in localized areas along the southeastern GrIS periphery 
and small values of less than 15 cm yr w.e. (units equivalent to g cm ' yr ) to the northeast 
and over the ice sheet plateau. The V 2 0 simulation also reflects the location of the equilibrium line 
along the southwestern side of the GrIS, local maxima northeast of the Uummannaq Fjord 
system and near Qaanaaq (Thule), and the sharply-defined ridge separating these two maxima. 
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Areas of ablation along the north and northeastern periphery are also represented, and both 
models indicate values less than 30 cm yr w.e. in a region extending southward from the 
plateau to near Tasiilaq. The GSN2 simulations indicate the large-scale features of GrIS SMB 
similar to RACM02 and ERA-I, including large values in the southeast, smaller values over the 
northeast plateau, and enhanced values inland of the ablation zone along the western side of the 
ice sheet. However the ablation areas along the western and northern sides of the ice sheet are 
larger than in either RACM02 or ERA-I, and more detailed features are largely absent. The 
effect of spatial resolution on the SMB distribution varies with location (Fig. 10). The 
southeastern maximum in SMB extends much farther into the interior of the ice sheet for the 
GSN2 field than for the higher resolution GSNA simulation. On the western side of the GrIS, the 
ablation zone is more narrowly defined and has greater magnitude at the higher spatial 
resolution, and there are larger values above the ablation zone. 

Averages computed over the GrIS as defined by the 1 .72x 10 km mask of 
Zwally et al. (2012) are shown in Table 1. The area encompasses GEOS-5 grid boxes of 
fractional land ice coverage along the GrIS periphery. For a particular GCM box along the 
margin, the fractional amount of area specified as land ice is dependent on the GCM grid 
spacing. Values obtained from only glaciated land surface tiles are indicated as “ice only”. The 
total value includes tiles within a GCM box that are not glaciated, and is used for direct 
comparison with other models having unknown fractional ice sheet coverage within a grid box. It 
may be seen that the average for the GSNV 2 simulation of 26.7 cm yr -1 w.e. compares favorably 
to the RACM02 estimate (27.4 cm yr -1 w.e.), while the 2° simulation value of 23.3 cm yr -1 w.e. 
is 15 percent less. For both model resolutions, the SMB computed using ice-only land surface 
tiles is smaller than the values computed from AGCM grid values. The ice-only values are 24.7 
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cm yr -1 w.e. for the GSN1/2 simulation and 18.2 cm yr -1 w.e. for the GSN2 integration. This is 
due to the larger amounts of runoff in peripheral ice-only tiles, as compared with AGCM grid 
values in which the runoff is averaged over both glaciated and non-glaciated land surface tiles. 
Table 1 shows a broad range of values for SMB components from RACM02, MAR, and ERA-I. 
For precipitation, the GEOS-5 simulations provide values of less than 40 cm yr _l w.e., which 
compare more favorably to ERA-I and MAR than to RACM02. The GEOS-5 integrations were 
averaged over the period 1980-2008, which is similar to the published averaging period for MAR 
and the same period as for ERA-I, while RACM02 was averaged over a 52 year period. The 
GEOS-5 simulations, the re-analysis, and the regional climate models are in rough agreement in 
finding evaporation has a small effect on the SMB. Interestingly, the GSN2 integrations find net 
deposition over the ice-only land surface tiles. As the peripheral Greenland topography extends 
over a wider area at coarser resolution, the northern interior regions of the GrIS are marginally 
more insulated from synoptic activity, and conditions are more conducive to sublimation. 

Finally, there is considerable uncertainty in the runoff term. The GSNA simulation value from 
the new surface scheme of 15.0±2.8 cm yr _l w.e. compares with estimates from RACM02 
(14.5 cm yr -1 ), MAR (10.1 cm yr -1 ) and ERA-I (17.7 cm yr -1 w.e.). 

A further evaluation of the SMB resulting from the new surface scheme is a comparison 
with the historical record of glaciologically-derived SMB estimates. An inventory of 331 
published values has been compiled by Bales et al. (2001), Bales et al. (2009), Cogley (2004), 
and van de Wal et al. (2012). A regression of these observations with values obtained from 
GEOS-5 integrations from interpolation of neighboring grid cells is shown in Fig. 11. Bales et al. 
(2001) separately listed historical records for years up to 1981 and more recent values. With 
notable exceptions, there is good agreement between the model SMB and the observations. Not 
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surprisingly, the GSNV 2 simulation compares more favorably at point locations than the GSN2 
run. The overall correlation is r = 0.563 for the GSN2 integrations and r = 0.648 for the GSN!4 
simulation. For the 2° model, there is a bias of -3.6 cm yr 1 w.e. with the observations, as 
compared to near zero bias for the V 2 0 simulation. But there is also a positive difference with 
observation at larger simulated SMB values from the GSN2 model as seen in Fig. 1 1(a). This is 
principally associated with the Bales et al. (2001) historical observations taken in the southern 
part of the ice sheet. 

Neither model resolution captures the large negative values tabulated in van de Wal et al. 
(2012). The observations were taken along the annually-visited K-transect over the period 1990- 
2010, at elevations ranging from 383 to 1850 m. At the lowest elevation site, SMB was 
consistently estimated to be less than -4 m yr -1 , while the modeled net ablation did not exceed 
1 m yr 1 for any of the K-transect locations. In a comparison with the CESM, Vizcaino et al. 
(2013) find similar differences in SMB with the K-transect. Vizcaino et al. (2013) attribute these 
differences to the location of these stations- in a tundra-surrounded fjord unresolved by the 
model- and to their location within the “dark zone”, an outcropping ice layer containing dust 
from an earlier period (Wientjes and Oerlemans, 2010). Excluding the observations of van de 
Wal et al. (2012), the correlation with in situ measurements is r = 0.671 for the GSN2 
integrations and r = 0.782 for the GSN'A model. 

A comparison of the simulated SMB has been made with ensemble members from ten 
similar AMIP models obtained from the fifth phase of the Coupled Model Intercomparison 
Project for the period 1980-2008 (CMIP5; Taylor et al., 2012). Not all of the participating 
CMIP5 models provided necessary variables- notably runoff- for computing SMB from the 
AMIP simulations. For some of the models, gaussian fields have been interpolated to an 
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appropriate regular grid. While the models shown here provide a runoff variable, it unclear from 
available documentation whether or not snow hydrology is explicitly simulated over the ice 
sheet, or if runoff results from peripheral non-glaciated land surfaces. The average SMB from 
models shown in Fig. 12 and Table 1 present a range of spatial resolutions and model skill. The 
coarse spatial resolution models of the Beijing Climate Center (BCC-CSM1.1) and the Canadian 
Centre for Climate Modeling and Analysis (CanAM4) depict only primary accumulation features 
of a maximum value along the southeast GrIS and a minimum region to the north. The 
remaining models excluding the those of the Geophysical Fluid Dynamics Laboratory (GFDL) 
approximate the 2° resolution of the coarse GEOS-5 simulations. Despite similarities in 
resolution among some of the models, the topographies are disparate, and are strongly related to 
both average SMB over the GrIS and local differences. Model topographies with a smaller area 
above 2000m (Table 1) such as the BCC-CSM1.1 have higher values for precipitation. None of 
the models with grid spacing larger than 1.3° capture the magnitude of local accumulation 
features suggested by RACM02. Higher resolution models from the GFDL present detailed 
patterns of coastal SMB, similar to the GSNA simulation and RACM02. In particular, the 
HIRAM-C360 SMB field indicates local maxima associated with topographic features along the 
western ice sheet margins. Both simulations indicate ablation regions in western coastal 
Greenland, but confined to the non-glaciated land surface. It may be seen from Fig. 12 and 
quantified in Table 1 that the RACM02 and ERA-I SMB fields are composed of ablation areas 
that are greater than 7 percent of the GrIS. Apart from the GEOS-5 models, only the Centre 
National de Recherches Meteorologiques (CNRM CM5) indicate similar ablation areas on the 
ice sheet. 
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The SMB for these AMIP model simulations, averaged for the GrIS (Table 1) range from 
25.0 to 46.5 cm yr -1 w.e. A large part of this range is attributable to differences in the amount of 
annual precipitation among the models, which ranges from 38.7 to 54.0 cm yr -1 w.e. But only 
three of the ten models examined indicate amounts of runoff greater than 5 cm yr -1 w.e. These 
are the CanAM4, the CNRM-CM5, and the MRI-CGCM3. The CNRM-CM5 utilizes a variation 
of the Douville et al. (1995) snow hydrology scheme over polar ice sheets, and attains values for 
runoff and SMB comparable to RACM02 and MAR regional models and the GEOS-5 
simulations. Interestingly, the Goddard Institute for Space Studies model (GISS-E2-R) provides 
a negligible amount of runoff over the GrIS but has the largest value among the AGCMs 
evaluated of 1 1.7 cm yr -1 w.e. for evaporation. In the averaged annual cycle of SMB (Fig. 13), 
summertime runoff is a dominant characteristic of the seasonal cycle of SMB in the re-analysis, 
and the GEOS-5 models, with each indicating negative values over summer months. The results 
are similar to those presented by Vernon et al. (2013) using regional models. Among the CMIP5 
models shown, only the CNRM-CM5 indicate negative SMB over summer months. 

Implications of differences in the representation of runoff on predictability are briefly 
summarized in Table 2, which gives the interannual trends in simulated SMB components. While 
the re-analysis is forced by satellite and in situ observations, interannual variability in AMIP 
models is primarily forced with SSTs and trace gas constituents (Gates et al., 1998). Along with 
the ERA-I, only the GEOS-5 simulations indicate a negative trend in SMB, which is marginally 
significant in GSN2 over for the integration period. In contrast the 10 models examined show 
significant positive trends in SMB, which is generally attributable to increases in the 
precipitation component. Precipitation trends in the AMIP models are larger than for ERA-I, 
which is larger than for the GEOS-5 models. The trend in runoff trend is significant and 
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dominates the SMB trend for the ERA-I re-analysis and the GEOS-5 models, while the runoff 
trend is as large in only a few of the other models. 

4. Summary and Discussion 

A new surface representation for glaciated land surfaces in the GEOS-5 model provides 

an improved near-surface temperature field and reduces biases in the net surface energy flux as 

compared to observations. In comparison to in situ AWS observations, a summertime air 

temperature bias -4.4°C is reduced to -1.4°C, and significant correction is also found in 

computed clear-sky surface temperatures in comparison to MODIS-derived values. The scheme 

eliminates a significant summertime bias in the net surface energy flux in comparison to in situ 

”2 

observations, and an erroneous annual mean net surface energy flux of up to 25 W m " found in 

the earlier method. These improvements reveal an underestimate of the surface longwave 

_2 

radiative flux of 4 to 6 W m that is compensated for by turbulent fluxes. A downwelling 
longwave bias is a chronic issue for models in polar regions (e.g., Ettema et ah, 2010a) and 
requires further evaluation, particularly with respect to cloud and aerosol properties. 

With the inclusion of a representation of snow hydrology and runoff, the primary 
variables for SMB are now computed by the GEOS-5 model. Spatial and temporal characteristics 
of the SMB simulated with a 14 ° grid spacing compare favorably with results using the 
RACM02 and MAR regional climate models, and with historic in situ observations derived from 
glaciological methods. Simulations of GEOS-5 at 2° grid spacing reproduce principal elements 
of the SMB field, but lack detailed features of accumulation along the western margins. In 
comparison to observations, the lower resolution simulations tend to overestimate the SMB in 
the southern GrIS by allowing maxima features in accumulation to extend farther inland, while 
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the average SMB over the whole ice sheet is less than that of the Vi° simulations. Other 
contemporary AGCMs examined at resolutions of 2° or coarser also have limited capability in 
simulating SMB beyond the primary features of large values along the southeastern coast and 
smaller values to the north and over the interior plateau. This suggests that a grid spacing of 
greater than 2° or about 200km is insufficient for adequately resolving coastal accumulation 
features. 

In comparison to the GEOS-5 model, an analysis of available CMIP5 AMIP simulations 
finds only a few models produce significant summertime runoff, a condition that marginally 
differs from that described by Randall et al. (2007). However, several recent studies have 
highlighted current work in this area. Using simulations of the GrIS with the Institut Pierre- 
Simon Laplace (IPSL) global model, Punge et al. (2012) describe the application of a multilevel 
snow hydrology model that includes compaction processes and prognostic albedo. Vizcaino et al. 
(2013) describe the surface configuration for 1° simulations of the CESM coupled model, which 
include snow hydrologic processes of compaction, percolation, refreezing and prognostic surface 
energy computations. Using a downscaling procedure, the CESM was found to compare closely 
to the RACM02 regional climate model. The SMB for the CESM plotted at 1° resolution is 
qualitatively similar to that of the GEOS-5 V 2 0 simulation in showing ablation areas along the 
western and northern edges of the ice sheet. 

Due to the presence of rough topography, spatial resolution is an important factor for 
adequate simulations of the GrIS. Nevertheless, the treatment of surface hydrological processes 
is also critical for reproducing the spatial variability and seasonal cycle of SMB. Additional 
issues of importance include the representation of spatially varying surface albedo, adequate 
vertical resolution to capture subsurface temperature gradients (Fig. 7), and use of an accurate 
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topography (e.g., Fig. 12; Box and Rinke, 2003). An adequate representation of the spatial and 
temporal characteristics of the surface temperature and SMB over the GrIS are seen as essential 
for coupling the AGCM with a dynamical ice sheet model. If potential contributions to future 
eustatic change from polar ice sheets are to be examined with GCMs, then current shortcomings 
in modeled SMB will need to be remedied. 

The surface representation described here will be used in a model coupling strategy to 
evaluate the relative importance of surface processes versus ice sheet dynamics in the evolution 
of the GrIS. The scheme is currently being implemented as part of numerical analyses and 
modeling conducted by the GMAO. Along with a parameterization for sea ice albedo, it will be 
incorporated as part of a group of improvements for polar processes in the upcoming MERRA2 
reanalysis. 
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List of Tables 


1 . SMB components precipitation (P), evaporation (E), and runoff (R) for the GrIS as defined by 
Zwally et al. (2012) from AMIP simulations of GEOS-5 and ten models of the CMIP5 project 
for the period 1980-2008, corresponding reanalysis values from ERA-I, and regional climate 
model output, in cm yr H w.e. The standard deviation of ensemble-averaged annual values is 
indicated in parentheses. The fractional ice sheet area above 2000m elevation from simulation 
topography, and the fraction area of ablation are given in percent. 

2. The trend in SMB components from AMIP simulations of GEOS-5, ten models of the CMIP5 
project for the period 1980-2008, and corresponding reanalysis values from ERA-I, in cm 

yr 1 w.e. The standard error of the ensemble-averaged trend is indicated in parentheses. 
Significant trends are shown in bold. 
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List of Figures 

1 . Schematic indicating the representation glaciated ice surfaces in the GEOS-5 model: (a) in 
original state and as represented in MERRA reanalysis, and (b) with application of 
snow model and ice representation. 

2. Greenland topography used for (a) 2.0° x 2.5° grid spacing and (b) V° x 5 / 8 ° grid. The contour 
interval is every 200m, with heavier line and labels for every 1000m contour . In (b), AWS 
are indicated with and “X” for GC-Net stations and “O” for DMI stations. Place names used 
in text are labeled. 

3. The net surface energy flux (a) and surface radiation (positive downwards) and turbulent heat 
flux components (positive upwards) (b) from GEOS-5 V 2 0 AMIP simulations in comparison to 
2000-2002 observations at Summit (72° 35'N, 38° 30'W; Hoch, 2005), and corresponding net 
flux (c), surface radiation, and turbulent heat flux components (d) for the S9 station within the 
ablation and percolation zone along the K-transect (67° 3'N, 48° 14'W; van den Broeke et al., 
2008). Model fluxes are averaged for the period 1980-2008. Bold solid lines indicate observed 
station values, dashed lines correspond to the CNTRLV 2 integration, and dotted lines indicate 
the GSNV 2 simulation, in W m 2 . 

4. The 1980-2008 averaged July surface albedo from GEOS-5 GSNV 2 simulation. The contour 
interval is 0.04. The margin of the ice sheet (Zwally et al., 2012) is indicated with a white 
line. 

5. Comparison of average Greenland AWS near-surface air temperatures from fifteen interior 
GC-Net locations (top) and twelve coastal DMI stations (bottom) with corresponding 
CNTRLV 2 and GSNA simulated values, in °C. Station locations are indicated in Fig. 2. 
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Model temperatures are taken from the nearest grid point and corrected to station elevation 
using the dry-adiabatic lapse rate. Diagonal lines denote exact correlation and ±3 °C. 

6. Average June- July-August clear-sky surface radiance temperature from GEOS-5 V 2 0 AMIP 
simulations minus MODIS temperature values for (a) CNTRLV 2 and (b) GSNA for 2000- 
2008, in °C. 

7. The sub-surface seasonal temperature cycle from (a) GC-Net AWS thermocouple temperature 
string at TUNU-N (78° 01'N, 33° 59'W), (b) simulation values for TUNU-N from GSNA, 

(c) GC-Net AWS values at Swiss Camp (69° 34'N, 49°20'W), and (d) simulation values for 
Swiss Camp from GSNA, in °C. Values are averaged for the period 1996-2008. 

8. Pan-Arctic difference of GSN'A AMIP simulation minus CNTRLA for June- July-August 
(a) 2-m air temperature, in °C and (b) sea level pressure, in hPa. Hatched areas indicate 
significance from the student’s t at the 95 percent level. 

9. Average surface mass balance from (a) RACM02 for the period 1958-2007, (b) ERA-I for the 
period 1980-2008, (c) GSNV 2 simulation for 1980-2008, and (d) GSN2 simulations for 1980- 
2008, in cm yr _1 w.e. The ice sheet margin is indicated with a red line. Topography is 
indicated with dashed contours for every 200m. 

10. The difference of GSN2 minus GSNV 2 average surface mass balance for 1980-2008, 
in cm yr 1 w.e. 

11. Comparison of 1980-2008 average SMB interpolated from neighboring grid cells from (a) 
GSN2 and (b) GSN'A simulations of GEOS-5 with corresponding observations from historical 
and recent values from Bales et al. (2001), PARCA values from Bales et al. (2009), Cogley 
(2004), and Wal et al. (2012), in cm yr -1 w.e. See text. The diagonal line denotes a 1 : 1 
comparison. 
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868 12. Annual mean Greenland surface mass balance from ensemble averages of CMIP5 20th 

869 Century AMIP simulations from 10 models, in cm yr H w.e. The average corresponds to the 

870 period 1980-2008. The ice sheet margin is indicated with a red line. Topography is signified 

87 1 with dashed contours for every 200m. 

872 13. Monthly (a) surface mass balance and (b) runoff for the GrIS from ensemble averages of 

873 CMIP5 20th Century AMIP simulations from GEOS-5 at two spatial resolutions, ERA-I 

874 reanalysis, and 10 models, in cm yr _1 w.e. The average corresponds to the period 1980-2008. 
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Table 1. SMB components precipitation (P), evaporation (E), and runoff (R) for the GrIS as 
defined by Zwally et al. (2012) from AMIP simulations of GEOS-5 and ten models of the 
CMIP5 project for the period 1980-2008, corresponding reanalysis values from ERA-I, and 
regional climate model output, in cm yr 1 w.e. The standard deviation of ensemble-averaged 
annual values is indicated in parentheses. The fractional ice sheet area above 2000m elevation 


from simulation topography, and the fraction area of ablation are given in percent. 



No. Ens. / 
Grid 
Spacing 

Area 

above 

2000m 

[%] 

P 

E 

R 

SMB 

Ablation 
Area [%] 

GEOS-5 

1 / Vi° 

61.1 






GSN'/ 2 (total) 



40.3(3.8) 

0. 8(0.2) 

12.8(2.6) 

26.7(4.4) 

10.6 

GSNA (ice only) 



39.8(3.7) 

0.1(0.1) 

15.0(2.8) 

24.7(4.5) 

9.6 

GEOS-5 

2/2° 

59.1 






GSN2 (total) 



34.5(2.5) 

1.6(1. 5) 

9. 8(1. 2) 

23.3(2.7) 

11.0 

GSN2 (ice only) 



33.7(2.5) 

-0.5(0. 1) 

16.0(1.9) 

18.2(3.3) 

12.1 

RACM02* 

-/1 1km 

61.0 

43.4 

1.5 

14.5 

27.4(2.4) 

13.8* 

MAR* 

- / 1 5km 


35.7 

0.2 

10.1 

25.3 


ERA-I 

-/ 0.75° 

58.0 

36.6(3.4) 

1. 4(0.2) 

17.7(2.7) 

17.6(3.8) 

7.9 

BCC-CSM1.1 

3/3.6° 

20.9 

54.0(2.8) 

5. 9(0. 4) 

3. 0(0. 2) 

45.2(2.7) 

0.0 

CanAM4 

4/3.6° 

41.7 

39.3(2.4) 

2. 1(0.2) 

12(0.9) 

30.0(1.8) 

0.0 

CNRM-CM5 

1/1.875° 

53.2 

38.7(5.3) 

2. 4(0. 2) 

11.2(2.1) 

25.0(4.6) 

7.7 

GFDL-HIRAM-C 1 80 

3/0.5° 

56.6 

43.7(2.6) 

2. 1(0.1) 

0.4(0. 0) 

41.2(2.5) 

0.0 

GFDL-HIRAM-C3 60 

2/0.25° 

57.5 

42.4(3.3) 

2.0(0. 1) 

0.2(0. 0) 

40.2(3.2) 

0.0 

GISS-E2-R 

6/2° 

55.7 

40.5(2.4) 

11.7(0.4) 

0.2(0. 1) 

28.7(2.1) 

2.5 

IPSL-CM5A-MR 

3/1.3° 

60.5 

41.7(3.2) 

3. 0(0.2) 

0.6(0. 1) 

38.1(3.1) 

0.0 

MPI-ESM-LR 

3/2.5° 

44.3 

42.7(2.8) 

0.2(0. 1) 

0.6(0. 1) 

41.9(2.7) 

0.0 

MRI-CGCM3 

3/1.5° 

53.9 

52.2(3.8) 

-0. 3(0.2) 

6. 0(0. 7) 

46.5(3.7) 

0.1 

NorESMl-M 

3/1.875° 

45.4 

52.1(2.8) 

2.1(02) 

1. 7(0.3) 

47.8(2.6) 

0.0 


*For 1957-2008, from Ettema et al. (2009). 

'Computed from interpolated 5km grid. 
^For 1990-2010, from Franco et al. (2012). 
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Table 2. The trend in SMB components from AMIP simulations of GEOS-5, ten models of the 
CMIP5 project for the period 1980-2008, and corresponding reanalysis values from ERA-I, in 
cm yr H w.e. The standard error of the ensemble-averaged trend is indicated in parentheses. 
Significant trends are shown in bold. 



No. Ens. 
/ Grid 
Spacing 

P 

E 

R 

SMB 

GEOS-5 

GSN'/ 2 (total) 
GSN’A (ice only) 

1 / V 2 0 

0.07(0.08) 

0.07(0.08) 

0.00(0.0) 

0.00(0.0) 

0 . 14 ( 0 . 05 ) 

0 . 16 ( 0 . 06 ) 

-0.07(0.10) 

-0.09(0.10) 

GEOS-5 

GSN2 (total) 
GSN2(ice only) 

2/2° 

0.01(0.06) 

0.00(0.06) 

0.00(0.00) 

0.00(0.00) 

0 . 09 ( 0 . 02 ) 

0 . 16 ( 0 . 03 ) 

- 0 . 08 ( 0 . 06 ) 

- 0 . 16 ( 0 . 07 ) 

ERA-I 

- / 0.75° 

0 . 13 ( 0 . 07 ) 

0.00(0.00) 

0 . 23 ( 0 . 04 ) 

- 0 . 09 ( 0 . 08 ) 

BCC-CSM1.1 

3/3.6° 

0 . 16 ( 0 . 06 ) 

0.00(0.01) 

0 . 01 ( 0 . 00 ) 

0 . 15 ( 0 . 05 ) 

CanAM4 

4/3.6° 

0 . 16 ( 0 . 05 ) 

-0.01(0.00) 

0 . 06 ( 0 . 02 ) 

0 . 11 ( 0 . 04 ) 

CNRM-CM5 

1/1.875° 

0 . 29 ( 0 . 11 ) 

0.00(0.00) 

0 . 13 ( 0 . 04 ) 

0 . 15 ( 0 . 10 ) 

GFDL-HIRAM-C 1 80 

3/0.5° 

0 . 17 ( 0 . 05 ) 

0 . 01 ( 0 . 00 ) 

0.00(0.00) 

0 . 16 ( 0 . 05 ) 

GFDL-HIRAM-C360 

2/0.25° 

0 . 23 ( 0 . 06 ) 

0.00(0.00) 

0.00(0.00) 

0 . 22 ( 0 . 06 ) 

GISS-E2-R 

6/2° 

0 . 17 ( 0 . 04 ) 

0 . 03 ( 0 . 01 ) 

0.00(0.00) 

0 . 14 ( 0 . 04 ) 

IPSL-CM5A-MR 

3/1.3° 

0 . 17 ( 0 . 06 ) 

0.00(0.00) 

0.00(0.00) 

0 . 16 ( 0 . 06 ) 

MPI-ESM-LR 

3/2.5° 

0 . 18 ( 0 . 05 ) 

0.00(0.00) 

0.00(0.00) 

0 . 18 ( 0 . 05 ) 

MRI-CGCM3 

3/1.5° 

0 . 19 ( 0 . 08 ) 

0.00(0.00) 

0 . 06 ( 0 . 01 ) 

0 . 13 ( 0 . 08 ) 

NorESMl-M 

3/1.875° 

0 . 20 ( 0 . 05 ) 

0 . 01 ( 0 . 00 ) 

0 . 02 ( 0 . 01 ) 

0 . 18 ( 0 . 05 ) 
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